Charge pair hopping and Bose-Einstein condensation in underdoped Mott insulators 
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Recently, a renormalized Hamiltonian has been derived by continuing spin states from the Mott 
limit of the t-J model to the underdoped region. We show that it naturally leads to a pairing 
mechanism in which the pair has a dual character. Its spin part is a spinon singlet which condenses 
at T* . The charge part is a real-space holon pair formed at T v < T* , which undergoes Bose-Einstein 
condensation at T c < T p . While neither is observable separately, the combination is. The mechanism 
is consistent with the small superfluid density, the decline of T c at small doping, and the existence 
of pairs above T c in cuprates, indicated by the observation of diamagnetism and Nernst effect. 

PACS numbers: 74.20.Mn, 74.72.-h, 71.27.+a, 74.20.-z. 



In a conventional superconductor the metallic state 
is characterized by quasielectrons, which form spatially 
overlapping Cooper pairs. In contrast, the elementary 
excitations in the metal phase of a hole-doped cupratc 
superconductor are hidden, and their nature unknown, 
which raises the questions: what is paired and how? Here 
we address this issue using the t-J model, which describes 
electrons of concentration 1 — x hopping on a lattice such 
that no site is doubly occupied. At half filling (x = 0), 
the system is a Mott insulator, with electrons localized 
as moments (or spins) which interact antiferromagneti- 
cally. Unexpectedly, the metallic state in cuprates is two 
dimensional (2d), even though these are (layered) 3d ma- 
terials. Anderson argued that in 2d electrons remain lo- 
calized even for x > 0, metallic conduction results from 
the motion of positively charged spinless holons. It is 
then natural to invoke continuity and construct a theory 
of the doped region by continuing the spin states from 
the insulator Experimental support comes from the 
fact that the carrier concentration is x, not 1 — x. Also, 
despite the reappearance of nodal quasielectrons below 
T c , the T = superfluid density p s (0) ~ x @, which is 
suggestive of holon pairing. 

In the t-J model, hopping by a localized electron can be 
described as an exchange of a spin- 1/2 particle (spinon) 
and a spinless holon, and the AF interaction as an ex- 
change of two spinons. This reflects the U(l) gauge sym- 
metry of the model. The electron field is then represented 
as c| = b\ a hi, where b icr creates a spinon of spin a and hi 
destroys a holon, such that the particle number at each 
site i is conserved and equals 1. The ground state at half 
filling is known: it is a mixture of the Arovas-Auerbach 
valence-bond (VB) state, in which bosonic spinons are 
paired into singlets, and the Neel state [H, Hj]. The VB 
state survives the destruction of AF order up to a tem- 
perature T*. However, despite many attempts, the con- 
nection with the doped phase (with fermionic holons) has 
not been established theoretically. 

Recently, we have solved the continuation problem and 
derived a renormalized Hamiltonian valid for small x , as- 
suming that hole motion prevents long-range magnetic 
order beyond some small x c 0, 0| • This implies that spin 
excitations are gapped since spinons are bosons. The 



new Hamiltonian has the symmetry of the original model, 
which allows us to determine its phases by continuing 
all the contiguous spin phases from half-filling. The ap- 
proach is blind: no experimental inputs are introduced 
by hand. Thus, we avoid the universal practice of as- 
suming two dimensionality of the normal state; instead 
it emerges as a consequence of the theory [f|. This is 
fundamental since 2d confinement is thought to be re- 
sponsible for the unusual behavior. It provides a serious 
test since it is almost impossible to confine metallic con- 
duction in the presence of 3d hopping; it has not been 
demonstrated in previous theories. 

Initially no new order is introduced, so that the phases 
are fully constrained by the symmetries at half filling. 
Yet, we obtain exactly two normal phases for small x as 
in cuprates Q, with properties that match experiments, 
but disagrees with other theories (see, @ for details). 
(1) Recent experiments have shown that, for small x, the 
high-T strange-metal phase actually behaves like an insu- 
lator, with no Drudc peak and a resistivity exceeding the 
Mott maximum, consistent with dynamically localized 
charge carriers @- In our case, since the high-T phase 
has the symmetry of the Hamiltonian, it is automatically 
an insulator, as holons are localized by gauge forces. (2) 
For T < T*, the pseudogap temperature, holons form 
a spinless Fermi liquid of concentration x, in agreement 
with recent experiments [8T-ITT1]. (3) The renormalized 
Hamiltonian has a pair hopping term, which automat- 
ically leads to d-wave superconductivity. (4) There is 
strong evidence for a spin gap Q which causes a down- 
turn in the paramagnetic susceptibility below T° > T*. 

Here we take a closer look at pairing. We show that 
ours is necessarily a strong-coupling Hamiltonian which, 
together with low density, induces holons to form real 
space pairs below T„. These undergo a Bose-Einstein 
condensation (BEC)[12| below T c < T p , which is quali- 
tatively consistent with the observed downturn of T r at 
small x, and diamagnetism [l3| and Nernst effect [bH 
above T c . A key prediction is a composite gauge-invariant 
"Cooper pair", which is observable above and below T c . 
It is a well-defined excitation of momentum q and en- 
ergy w(q), whose charge is carried by a mobile (bound) 
holon pair, which picks up the local phase of spinon sin- 
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glcts (representing localized electrons) already condensed 
below T*. 

The t-J Hamiltonian for a layered 3d system is given 

by 

H = -5>i4<> - 2^./,...\;,A,.. (1) 

ij 

The first term describes electron hopping from j to i such 
that no site is doubly occupied, the second is the ex- 
change interaction between spins. Here, Ay = \\b^bj^ — 
bi]bj^] destroys a spinon singlet, and Jy = Atfj/U, 
where U is the Hubbard repulsion. Minimally, we need 
nearest-neighbor (t) and next-nearest-neighbor hopping 
(£') within each plane, and nearest-neighbor out-of-planc 
hopping t z << t. For cuprates, t/J ~ 3 — 4. In general, 
t'/f is small, so is J' = (t' /t) 2 J, and will be neglected. 

As described in Q, renormalization proceeds in two 
steps. Hole motion is violently opposed by AF correla- 
tions, which localizes the hole in a small region, within 
which it hops rapidly. The renormalized hole - bare hole 
plus surrounding spins - moves slowly through the lattice, 
preventing AF order. By continuity, the Hamiltonian has 
the same form as the original one, with renormalized pa- 
rameters t e ff, J e ff etc, except that spinons have a gap 
A s . The spin gap is the singlet breaking scale ~ 2A S . 
At small x, we expect J e //j to scale with J. One-hole 
calculations suggest that for t > J, in-plane hopping am- 
plitude t e ff < J, and also scales with J. The spin gap 
allows us to decouple unpaired spinons from hopping, 
and derive a minimal Hamiltonian involving sublattice- 
preserving hopping of renormalized holons and singlets 
only. Now, due to the gap, singlets form within the plane 
since J z /J = (t z /t) 2 « 1 and a spinon belongs to one 
singlet at a time. As a renormalized hole hops it breaks 
a singlet. For hops within the plane, the excess energy f2 
is removed if the hole makes a second hop, and the sin- 
glet is reconstructed. Eliminating the intermediate state 
perturbatively we obtain [(| 

H 1 = - t ^(l-x)J2^AiMhl, (2) 

ijl 

which describes nearest- and next-nearest-neighbor hop- 
ping within the same sublattice of the plane, accompanied 
by a singlet backflow. The (1 — x) factor reflects the ab- 
sence of a hole in the intermediate state. For a hop to 
a neighboring plane, two unpaired spins are created, one 
in each plane. Now the system can not relax by the hole 
hopping to a third plane since that would create more 
unpaired spins, which proliferate exponentially, blocking 
hole motion in the z direction. 

The system can also relax if a second hole follows the 
first, and the singlet is reconstructed. The net effect is 
hopping by a holon pair to a parallel link, accompanied 
by a singlet backflow. This mechanism also works for 
interplane hopping, since the entire singlet is transferred. 



Eliminating the intermediate state, we get 

ijlm;z ij;z 

(3) 

where t sz = t 2 zeff /4fl, and Cy = (c^c^ - c^Ci^)/2 = 

— AijFT - destroys a singlet made from physical electrons, 
which is equivalent to destroying a spinon singlet and 
creating a holon pair with F^ = hjhj. The first term de- 
scribes intraplane, and the second, interplane pair hop- 
ping. The full Hamiltonian is gauge invariant and has 
the symmetries of the original model, plus an additional 
symmetry: total number of holes in each sublattice is 
conserved. 

Since the high-T phase has full symmetry (nonorderd) 
holons are localized by gauge forces. A pseudogap metal 
appears below T* as the spinon singlets condense, allow- 
ing holons to propagate coherently. It is connected to 
the VB state at x = 0, and is characterized by the order 
parameter 

An =< Ay >= Ae l i Q -( ri - r i\ (4) 

or, its gauge-related copies, where Q is the two-sublattice 
wave vector. This state has a quantum lattice order: on 
average singlets connect spinons on opposite sublatticcs. 
If we replace Ay by its average Ay , the single- hole term 
[Eq. 2] describes free holon hopping within the sublat- 
tice, with a spectrum eh = —2th + 2th(sink x + sinfc y ) 2 + 
4D it s A 2 sin k x sin k y , where th = t s A 2 (\ — x). The last 
term is the Hartree contribution from the pair hopping 
term, and D\ < x is the average hole hopping amplitude. 
At low-T, the pseudogap phase is thus a spinless Fermi 
liquid of concentration x. The corresponding small Fermi 
pockets are not gauge invariant, but have been seen in- 
directly in de Hass-van Alphen type experiments fllj ] . 
The condensate part of the pair-hopping term becomes 

H 2h = -t s A 2 [ F U z ) F mi(z) + h-c] 

ij,lm;z 

-t zs A 2 Y,[FU^(z + l) + h.c\. (5) 

It clearly leads to 3d superconductivity via pair conden- 
densation, so that F ?J =< Fy >^ 0. The electron pair 
amplitude is then < C v] >= - < Ay >< F, y >^ 0. 
A T = MF analysis showed that it is a robust d-wave, 
essentially due to the symmetry of the VB state at x = 0. 

The MF approximation clearly would not work for 
T > since the intraplane pair-hopping energy scale 
t s is essentially the same as that for single-hole hop- 
ping, (t s (l — x)) (Eqs. 2, 3) because they arise from 
the same singlet breaking mechanism. It would induce 
holons to form real-space pairs, particularly since (a) at 
small x, pairs would not overlap, and (b) they can re- 
duce energy further by delocalizing in the z-direction. 
Here we show this by analyzing pair fluctuations using 
functional integral methods [l5j and focusing on T >T C . 
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A holon pair field is written as a two-component vec- 
tor field: Fij{r) = Fj^(r), where r\ = (x,y). Let 
F "' (r) = ^=E p e 4k -( ri +"/ 2 )- i "^, ) (k, W ), where V is 

the unit vector, and p = (k, tu). The Hamiltonian density 
is given by 

H = J2 mKh p - J2 E n W (p)Fv (p), (6) 
p pi; 



where £(k) = e h (k) - /U/j, and 

-E^.y (k) = 2t cos fc^ + 2t z0 cos fc z 



(7) 



is the hopping energy for a holon pair. Here to = 
t s j4 2 , and f z0 = t zs A 2 . We introduce an order- 
parameter field A v (p), and do a Hubbard-Stratonovich 
transformation to obtain the action S — So ~ 
E M [A;(P)A,W - v^„(k)A; (k)F v (p) + c.c.l where S 
is the free (quadratic) holon part. Integrating out the 
holons and keeping terms to second order in A we obtain 
the effective action describing pairing fluctuations: 

Sn = -5>„„ - (^(k) J B p (k)) 1 / 2 n w (p)]A;(p)A p (p), 

pr/p 

(8) 

where U vp (p) =< F*(p)F n (p) > is the pair correlation 
function for noninter acting holons, which is given by 



n w (q,w) 



1 E 



sin k v sin k p 



JV^iw-£(k + q/2)-£(k-q/2) 

[tanh e(k+c^ + tanh e(k-q/2) 1 



2T 



2T 



(9) 



The sum is over the 2d Brillouin zone and fcs = 1. 

The thermodynamic potential is given by = fix +^2, 
where tti is the free holon part, and ^ = T^ y ln[(l — 
A+(p))(l — A_(p))], is the fluctuation contribution, with 



' Fy^yy {{Fx^xx Fy^yy) 



~ 4 X Ey II ^ y II y X } 



(10) 



Since E x (0) = E y (Q) = 2t + 2t z0 , we obtain, for p = 
(noting that U xy < 0) A± = £ K (0) [11^(0) T n xy (0)]. 
The order-parameter equation is given by 1 = A+(0), the 
larger eigenvalue. At T c it reads 



1 



to + t z Q 



This corresponds to a d wave. The holon density is given 
by x = x\ + X2, where 



i = ^E[ 1 - tanh (e(k)/2T)] 



2N 

k 

is the contribution from free holons, and 
x 2 = -N- 1 da 2 /dfi h 



(12) 



(13) 



is the fluctuation contribution. 

Solving Eqs. (11-13) we obtain /i^ and T c for small x. 
We take the bottom of the holon band to be at zero. In 
the MF (BCS) approximation fluctuations are neglected 
(x = xi). The MF T c , which we denote by T p , is a 
measure of the pair binding energy, and remains finite as 
x — > 0, as shown in in Fig. (1). When fluctuations are in- 
cluded, real-space pair states appear as poles of the two- 
holon Green's function, i.e., as zeroes of Re(l — A±(q, u>), 
which arc at uj = oj±(q) = E p ,±(q) — /j, p . (We have made 
an analytic continuation to real frequencies). Here, E Pt ± 
is the energy of the pair, fi p is the pair chemical poten- 
tial, such that w±(q) > 0. We need to consider only the 
lower pair band E p+ (q) which has states below the free 
holon-holon {h-h) continuum, therefore are long lived. 
Once formed, the states near the bottom of the band 
(which is at q = 0), arc quickly occupied by holons so 
that X\ is negligible, and x sa x%- Then, ~ E p+ (0)/2 
is negative, and binding energy is ~ 2|/i/j|. Bose con- 




FIG. 1: Critical temperature T c and the holon pair dissoci- 
ation temperature T p as a function of hole density. In the 
region, T p > T > T c , pairs coexist with holons. 

densation occurs for ui = w + (q = 0) = 0, so that 
Hp = E p+ (0) ~ 2\fj,h\- Numerical solution shows that, 
for small x, \fih\ ~ t a » T c , so that pairs are strongly 
bound. Since pairs are bosons we need to consider only 
states near q = 0. It is easily shown that X± is symmetric 
in the q x , q y plane. Then for small q the pair spectrum 
has the form 



which gives 



w+(q) = a(q 2 x + ql) + bql 



T c a 4.66(z 2 ) 2 / 3 (a 2 &) 1/3 - 



The parameters a, b also depend on T c , and indirectly on 
x, and obviously x\ — x — xi is not zero. However, these 
corrections arise from free holons, which require a finite 
energy = 2|/i/j| to produce. For small T c /\fXh\, we obtain 
xi w er\^l T 'Zx, where Z x = ± J2 k e~ eh{W)/Tc ■ We find 
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from numerical calculations that a ~ to, and b scales is 
tzo << to- Since, T c also depends on X2 ~ x, T c j p^ ~ 
(xH^/to) 1 ^ 3 is small, and aq vanishes exponentially. Fig. 
1 shows T c and T p function of x. 

The bound pairs continue to exist up to T p > T c . 
Holons are gapped below T p , but there is no sharp tran- 
sition. Now, a mobile holon pair is hidden, as it is not 
gauge invariant. However, in the present case, it becomes 
observable through the gauge-invariant physical electron 
(singlet) pair represented by: Cy = —Ay Ft.. Ordi- 
narily the electron pair Green's function G e ^~, vair (p) = 
— < C v (p)C*, (p) > is incoherent since it is a convolu- 
tion. However, a coherent part emerges as the spinon 
singlets condense below T* . Replacing Ay by its mean- 
field value, we obtain G^J pmr {p) = A 2 G* „,(—p), where 
G vv>(p) = ~ < F r,{p)F*'(P) > To determine G F we use 
the exact relation (E^E^) 1 / 2 < F*F n , >= —5 m i+ < 
A*A V , >. This leads to 

G e J:Z pmr (p) = — — ^ nn -. + incoherent part, (14) 
" iuj + uj + (— k) 

which is a key prediction. Here 

W {E v E v ,y/2(d\ + /duj) 

evaluated at ui — w(k); u±(rf) are the eigenvectors 
which diagonalize the effective action [Eq. 8], yielding 
cigncvalues A±. Eq. (14) describes a well-defined physi- 
cal (gauge-invariant) pair excitation. However, unlike a 
Cooper pair, it is observable above and below T c , and 
repesents both electrons and holes. Since the pole is 
at negative energy — w+(— k), the actual excitation rep- 
resents a mobile hole pair of momentum k and energy 
w_l_(k), of charge 2e, The spin part consists of a spinon 



singlet (representing a localized electron pair). These are 
already condensed below T* , and tend to overlap because 
of high density (1 — x). With increasing x and/or T, the 
excitation will be less sharp due to broadening by phase 
(gauge) fluctuations associated with Aij, and the factor 
A 2 , which vanishes for T >T*. 

These results resolve a number of puzzling issues. 
Holon pairing implies that p s (0) ~ x. Superconductiv- 
ity by BEC explains the decline of T c with decreasing 
x. That T c does not vanish at some xsc ~ x c is not 
suprising since the renormalized Hamiltonian is not valid 
near x c , where the spin gap is small, AF correlations are 
longer ranged and compete with superconductivity. Ex- 
perimentally, T c cx pi, i.e., both vanish at same xs c , as 
if the effective carrier density is x — xsc- Early reports 
indicated that 7 = 1 1 1611 : however, recently 7 = 0.61 
has also been reported [l7[, which is not far form 2/3. 
The existence of charged bound pairs above T c has been 
shown to account [l|| for the observed diamagnetism [l]| . 
On the other hand, since the spin part of the pair is con- 
densed below T* > T c , and the associated phase (gauge) 
fluctuations couple to holon pairs, vortex type excita- 
tions are likely to exist. These may account for the ob- 
served Nernst effect [l4[, although more work would be 
needed to sort out these issues. Interestingly, inthis the- 
ory, spinon singlet condensation is two-dimensional, so 
that the transition at T* is probably Kosterlitz-Thouless 
type, hence the lack of singularities. Then phase fluctua- 
tions are also two-dimensional. The theory is consistent 
with the observed decoupling between spin and charge 
responses. As T decreases, the paramagnetic suscepti- 
bility decrease starting at the spin-gap temperature T , 
and seemingly unaffected by charge pairing below T p and 
superconductivity at T c [1, Q • Finally, though not con- 
sidered here, nodal electrons can appear as collective ex- 
citations [19( , but would not directly take part in pairing. 
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